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Abstract 

The full problem of flame stability for the two-reactant model, which 
takes into account thermal expansion effects for all disturbance wave lengths, 
is examined. It is found that the stability problem for the class of two- 
reactant flames is equivalent to the stability problem for the class of one- 
reactant flames with an appropriate interpretation of Lewis numbers. 
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1. Introduction. Recently, Jackson and Kapila (1984, 1986a, b) reexamined the 
problem of the stability of a plane, premixed flame in the near- 
equidiffusional (NEF) framework. The study was aimed at understanding the 
role of thermal expansion on flame behavior. Neutral stability curves were 
obtained in the Lewis number-wave number plane and contrasted with the 
classical results of Sivashinsky (1977), who had employed the Constant Density 
Approximation (CDA). 

The analysis was based on the simple kinetic model A B which, as is 
well known, is also appropriate for a multi-reactant system when a single, 
deficient component controls the combustion process. Near stoichiometry, 
however, revisions are necessary since the number of controlling reactants is 
more than one. The first effort in this direction was due to Sivashinsky 
(1980) who examined the stability of a two-reactant model at EXACT 
stoichiometry, under the CDA; Joulin and Mitani (1981) later extended the 
result to near-stoichiometric mixtures. The essential point made by these 
studies is that the single-reactant results hold even near stoichiometry, 
provided the Lewis number is replaced by a suitably defined effective value. 

The purpose of this note is to show that even when the CDA is abandoned 
and thermal expansion fully taken into account, the notion of an effective 
Lewis number persists. Thus, upon suitable interpretation, the results of 
Jackson and Kapila (1984, 1986a, b) can be applied to two-reactant systems. 


2. Governing Equations. Consider the two-reactant model 


V l^l + V 2^2 + Products, 
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where and M 2 denotes the chemical symbols of the deficient and abundant 
reactants respectively. Here is the stoichiometric coefficient for 
species j . We take the reaction to be irreversible and of Arrhenius-type so 
that its rate expression is of the form 

(2.1) u = B p 2 Y^ 11 ex P (-E/RT), 


where p is the density, Yj and Y 2 are the mass fractions of the deficient and 
abundant components, respectively, and T the temperature. The parameters 
appearing in the rate expression are the activation energy E, the universal 
gas constant R, the pre-exponential factor B, and the reaction orders m and n 
with respect to the deficient and abundant reactants, respectively. 

The basic equations of the two-reactant system can be modelled by the 
following equations: 


( 2 . 2 ) 


a. + p7 . v -o , p-sjr. 


(2.3) 


Dv 

p df + vp = y{ i V(7,v) + v2y} » 


/V /V 


(2.4) 


r DT DP 
P p Dt “ Dt 


= X V^T + Qw, 


(2.5) 


DY j 

VS 


VjWjU) , j = 1,2. 


Besides Yj and T, the state of the system is determined by the density p , 
pressure P and velocity y . The other quantities appearing above are the 
specific heat at constant pressure Cp, viscosity y , thermal conductivity 
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X , species diffusion coefficient Dj of species j, chemical heat release Q, 
molecular weight Wj of species j and molecular weight W of the inert. 
Finally, D/DT denotes the material derivative, defined by 


( 2 . 6 ) 


D 

Dt 


3_ 

3t 


Z- = + v7 


Equations (2.1) - (2.5) are supplemental by appropriate boundary conditions 

and jump conditions across the reaction zone. 

For nondimensionalization of the governing equations, the state of the 

fresh mixture is chosen as the reference, defined by the constant values 

P , p , T and v , while the initial concentrations are rendered 
o o o o 

dimensionless relative to the initial concentration of the deficient reactant 
Y^q. The characteristic length and time scales are taken to be the adiabatic 
diffusion scales 


(2.7) * 0 " X/<p o C p v o ) ’ t o = X/(p o C p v o 2) * 

It is convenient to introduce the flame Mach number as M q = V Q / C 0 > where c 0 

is the frozen sound speed in the fresh mixture. Finally, pressure deviations 

2 

from the ambient are referred to yM Q P , y being the specific-heats ratio. 
The governing equations then take the dimensionless form 

( 2 . 8 ) 

(2.9) 


jjj: + P V»V = 0 , 1 = pT , 

/V 

Dv 

p dT + vp = Pr (j V(v»v) + v 2 v} , 


P 


DT 

Dt 


( 2 . 10 ) 


V 2 T + crtl , 
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( 2 . 11 ) 


DY j 

p dT" = Lj 1 v 2 Yj - , j - 1,2 , 


( 2 . 12 ) 


= D v up Y^*" 1 p 2 Y“Y%- 0/T . 


M 


2 1 To 10 


12 


The yet undefined nondimensional parameters appearing above are: 


Pr = y C /X 
P 

“ - QY 10 /C pVl“l 

0 = E/RT 

o 

e = (1 + o) 2 /0 


Lj - X/p Dj c p 

D = AB/p C c 2 
o p o 

a . = v .W. /v , W. 
J 3 3 11 


Prandtl number, 

Heat release number, 

Activation energy, 

Reciprocal activation energy, 

Lewis number for species j , 
Damkohler number, 

Parameter involving stoichiometry. 


When a steady, planar, adiabatic flame is analyzed in the limit 
0 ■* «, i.e., e + 0, the following expression for the Mach number M q is 
obtained (see, e.g., Sen and Ludford (1979) and Mitani (1980)): 


(2.13) 


where 


D 



o 


q m+n+1 (l-ta) 2 
2 G(m,n,A) 


-1-m-n 1-m-n 

e Y 10 

TT T 111 T 11 

V l W l P o L l L 2 


0 

1+a 

e (1 + 0( e)} , 


(2.14) 


G(m,n,A) = /“ £ m (5+A) n e -5 d? . 


Here, A is the normalized final concentration of the abundant, given by 


(2.15) A = (4.-1) 

and <f> is the equivalence ratio ^ v i^i^ 20^ V 2^2^10^ * Equation (2.13) is valid 
for the entire range of mixture ratio. However, we are only interested in the 
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mixture-ratio range $ - 1 = 0( e ); i.e., near-stoichiometric mixtures, since 

the reaction rates are then governed by the deficient component as well as the 

abundant component of the mixture. The special cases of very-near 

2 

stoichiometry (<fr - 1 = 0( e ); i.e., A 0) and away-from stoichiometry 

(<J> - 1 = 0(1); i.e., A ») can be considered as limiting cases of the near- 

stoichiometric analysis, and will be discussed in the conclusions. Equation 

(2.13) also shows that the Mach number is exponentially small in the limit of 
large activation energy. 

The near-equidif fusional formulation (NEF) , in addition to large 
activation energy, is characterized by near-unity Lewis numbers and nearly- 
uniform enthalpies, i.e., 

—1 ^1 

(2.16) L J = 1 - e (^i) , T + aY x = H = 1 + a + eh, 

(2.17) L-* - 1 - «£ > , T + 2- Y 2 E G - 1 ♦ Y 2£ ♦ eg. 

Here, = ^20^10 ^ enotes t ^ ie dimensionless mass fraction of in the 

fresh mixture. Then, to leading order in e , equations (2.8), (2.9) remains 
unaltered, while equations (2.10) and (2.11) are replaced by 

DT 2 

(2.18a) p— - V T = 0 in the burnt region, 

(2.18b) T = 1 + a in the unburnt region, 

(2.19a) 

(2.19b) 


p ^ = 7 h + <4> v ' T * 


P || - v 2g + (-£) v 2 t. 
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Gquations (2.8), (2.9), (2.18) and (2.19) now constitute the basic set. 

The following analysis is by now standard, and the reader is referred to 
Jackson and Kapila (1984) for details, hereafter referred to as I. For a two- 
dimensional flow, let the reaction-zone location be given by 
x = x^ = F(y,t). It is convenient to shift to a reference frame in which 
the reaction zone is stationary. Let a new coordinate C be defined by the 
transformation 

(2.20) x = F(y,t) + K , 

and let the unburnt and burnt regions be confined, respectively, to 
5 < 0 and 5 > 0. Then one obtains equations (8)-(22) of I, except now the 
enthalpy equation (141) is replaced by the enthalpy equations for h and g, 

£ 

(2.21a) p (h t + Sh^ + vh y ) = V 2 h + (-i-) V 2 T, 

£ 

(2.21b) p(g t + Sg^ + vg y ) = V 2 g + (-|) V 2 T, 

while the boundary conditions (211) and (221) at 5=0 are replaced, 

respectively, by 

h / 2 w 1 / rt 

(2.22) 6T = ^ {/ “ (A+g-h+5) n e" 5 d5} 1/2 , 

5 /(1+Fp G(m,n, A) 

(2.23) 61^ = - (^-) 6T ? , 6g 5 = - (-^) ST^. 

Equations (8I)-(22I), together with (2.21)-(2.23) describe completely the NEF 


problem 



3. Linearization. Small perturbations of the steady flame are governed 


by linearized equations, which can be obtained by setting 


(3.1) <J> » ij> 8 (5) + ot|i (5,y,t) 

for all variables and taking the limit o 0. We focus our attention only on 
the enthalpy equations (2.21) and the boundary conditions (2.22) and (2.23) 
since linearization of the other equations yields the identical equations of 
I. 

Linearization of (2.21) leads to 


(3.2a) 


p S h. + h + m hf = h + h - F hf 

t 5 5 55 yy yy 5 


* i 

+ (— )(T__ + T -FT®), 

a 55 yy yy 5 


(3.2b) 


P S 8 t + g 5 + m g 5 ■ hi + 8 yy - Vc 


+ (— ) (T__ + T -FT?) 

a 55 yy yy 5 


for h and g, respectively, while linearization of the jump conditions (2.22) 
and (2.23) at 5=0 yields 


(3.3) 6T 5 = " I {Tr+ H (g " h »> 

£ £ 

(3.4) 6hg = - (-^-) 6T 5 , 6 g ? - - (-i) 6 T^ . 

Here, 
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(3.5) 


H = n 


G(m,n-1 ,A) 
G(m,n,A) 


If we define a new function f by 


(3.6) 


f - h + H (g - h). 


then (3.2) reduces to a single equation for f, namely 


(3.7) 


p 8 ?. + f + mf® = f + f - F f ® 

t 5 5 55 yy yy 5 


+ - (A. + H (Z n - £.)) {T__ + T - F T® } , 
a ' 1 2 1 ££ yy yy K 


while (3.3) and (3.4) determines the jump conditions at 5=0, 


(3.8) 


6 T - - | f , 6 f g = " 4 (£ 1 + H ( Jl 2" £ l )) 6T ? 


These equations are identical to the linearized enthalpy equation for the 
premixed flame (see (301), (361), and (421)) with a reduced Lewis number 
defined by 


(3.9) 


l c = Aj + H U 2 - 4 X ). 


4. Discussion and Conclusions. The stability problem of the two-reactant 
flame consists of the linear system (24I)-(42I), with (301) and (361) replaced 
by (3.7), and the boundary condition (421) replaced by (3.8). In particular, 
we find that this system is governed only by the parameter Jt^, 


defined in 
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(3.9). It is only here that the mixture ratio appears through H, and thus we 
consider the following three limiting cases. For near-stoichiometric 
mixtures, H is given by (3.5) and can not be simplified further. For very- 
near stoichiometry, H is given by 


(4.1a) 


H = mTn + 0(e) ■ ♦ “ 1 - 0 (e 2 ), 


and the corresponding reduced Lewis number becomes 


(4.1b) 


mi. + n £» 

£ —7— + 0 (e), 

c m + n 


i.e., the average value of the two Lewis numbers weighted by their respective 
reaction orders (Joulin and Mitani (1981)). Finally, away from stoichiometry, 
H is given by 


(4.2a) H =^~ y + 0 (e 2 ) , *-1-0(1), 

with corresponding 


(4.2b) 


l m t* +0 (e ) . 
c 1 


Thus, for far-from stoichiometric mixtures, the stability problem reduces to 
that of a one-reactant flame governed by the Lewis number of the deficient 


component 
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Sivashinsky (1980) and Joulin and Mitani (1981) have shown that within 

the diffusional-thermal model, the stability problem for two-reactant flames 

at or near stoichiometry reduces to the stability problem for single-reactant 

flames with a reduced Lewis number* We have shown that the above result also 

holds for the more important case when thermal expansion is accounted for, and 

for the entire range of mixture ratio. The following general statement can 

now be made: The linear stability problem for the class of two-reactant 

flames is equivalent to the stability problem for the class of one-reactant 

flames with Lewis number L = 1 + e(l /a). The effects of heat loss and 

c 

gravity do not alter this conclusion. 

Acknowledgements. The author would like to thank A. Kapila for helpful 
comments in connection with this work. 
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